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Abstract. For a Hamiltonian system x = JH'(x) in R 2n its two-system is defined in the 
phasespace R 2n x sp(2ra, R)$ as follows: x = J (If' (x)+tr ($#"))', 6 = [$, H"(x) J]. In 
a sense, it is a combination of the original system and its system in variations with feedback. 
We study the Hamiltonian forms of the two-system and its analogs. In particular, we show 
a relation between the Poisson-Lie bracket on the dual space to Lie algebra sp(2n, R) and 
the canonical bracket. 



Introduction 

In this paper we describe one special construction in the theory of Hamiltonian system called 
the two-system. Our point of view is different from that in |BBK 1994j . where such systems 
were defined, and in [BK 1994J, where the term "two-system" was introduced (in the context 
of semiclassical approximation in quantum mechanics). Here we treat the two-system as a 
certain union of the given Hamiltonian system and its system in variations. Versions of the 
construction, their Hamiltonian structures, and general properties are discussed. 

As an example, the following fifth order system of ordinary differential equations with 
parameter e 

q = p, p = —(q + eq 3 ) - 3sqa, 

a = 2(3, $ = -(1 + 3eg 2 )a + 7, j = -2p(l + 3eq 2 ). 

is the two-system corresponding to the underlying system with one degree of freedom and 
Hamilton's function 

H{q,p) = — - — + e—. 
This system with e 1 was analyzed in jS 1994j by the normal form method. 

1. A hamiltonian view on the system in variations 

a. Consider a Hamiltonian system with n degrees of freedom 

i=JH'(x). (1) 

Here x = (x 1 , . . . ,X2 n Y denotes a point in the phase space R 2n as well as a vector-function 
x(£); H(x) is Hamilton's function, i?'(x) its gradient 
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J is the standard symplectic matrix: if I n denotes the unity matrix of order n, then 




b. Along with system ([TJ consider its system in variations 

y = J#"(x(t))y, y=(y 1 ,...,y 2n ) t , (2) 
where x(t) is a solution of system (£[]), and if"(x) is the Hesse matrix, 

H" = 

13 dxidxj 

The following is a simple and well-known fact. 

Proposition 1. Let x(t) be any function (not necessarily a solution of system (£[])). Then 
system ^ is a time- dependent Hamiltonian system with respect to variables y, with quadratic 
Hamilton's function 

H*(t), y) = iy^"(x(t))y. (3) 

c. Our intent is to combine systems and (J2J) into an autonomous Hamiltonian system. 
A straightforward union of and © suits to this purpose only in exceptional cases. 

Proposition 2. Suppose the function jF(x, y) is such that system has the form 

x = 4 

<7X 

and at the same time system (J2J) has the form 

T/ien -ff(x) zs a polynomial of degree 2 or /ess. 
Proof. It follows from and (0J) that 

ox 

Hence ^(x, y) = H (x) + C(y) with some function C(y) independent of x. Then (j^J) becomes 
y = JC'(y). It is a system independent of x, which doesn't agree with (J2J unless if" is a 
constant matrix. □ 

Proposition 2 implies that in general any Hamiltonian combination of systems ((IJ and 
(J2J) has to be artificial in the sense that it will alter the dynamics. Nevertheless such a 
combination may arise naturally from another perspective. The Hamiltonian combination 
we deal with in the paper, or rather its generalization called the two-system, first appeared 
in the context of semiclassical approximation in quantum mechanics [BBK 1994J. We do not 
discuss quantum mechanical applications here. 



2 



2. Hamiltonian combination in the vector form 

Construction. Consider the phase space IR 4n with coordinates (x, y) and a symplectic struc- 
ture given by the matrix J © J, so that (xi,x i+n ) and (yi,yi+ n ) are pairs of canonically 
conjugate variables. Define Hamilton's function in M 4n as the sum of Hamilton's functions 
of systems (JTJ and (0) 

W(x,y) = if(x) + iy*ff"(x)y. (6) 
Setting in (HJ), (0) J 7 = H, we obtain the Hamiltonian system 

* = j(in*)+~y'#''(x)y)\ (7) 
y = J#"(x)y- 

The prime (') denotes differentiation with respect to x throughout. 

Definition 1. We call system (J7|) the vector form (for a Hamiltonian combination of system 
and its system in variations (J2J) ))■ 



3. Hamiltonian combination in the matrix form 

a. Introduce the 2n x 2n matrix M = yy* and rewrite © in the form 

W(x, y) = i?(x) + ^tr(tf"(x)M). (8) 

The second equation in (JJJ) is equivalent to 

M = JiJ"(x) M + MH"(x) J 1 . (9) 

Set 

$ = JM, A = H"J. 
Using the identities J* = — J, J 2 = — i2n, write Q in the Lax pair form 

$ = [A, $] = A$ - $A 

The entire system becomes 

x = J (#(x) - ^tr(A(x) <&))', 
$ = [A(x), $]. 

Definition 2. System ^TDJ) with a vector-function x(t) and a matrix-valued function is 
called the matrix form (for a Hamiltonian combination of the systems (0) and (J2J ))• 

Remark. The vector form (j2J) is contained in the matrix form (jlOjl and corresponds to the 
invariant set of such (x, $) that M = — J$ is a symmetric non-negative matrix of rank 1. 
Section elaborates on this correspondence. 
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b. An arbitrary matrix M can be decomposed into its symmetric and antisymmetric parts 



M + M t M — M t 
M s = ~ 2 , M a = . (11) 

Now take M = — J$. Let us show that system (jlUj) respects the decomposition (fTTj) . and 
that only the symmetric part M s has an effect on the evolution of x. 

Proposition 3. Let (x(t), JM(t) ) be a solution of system iflUjl. M = M s + M a be the 
decomposition (fTTjl. Denote $ s = JM S , $ a = JM a . Then 

(i) (x(t), $ s (i) ) is a solution of system (fTU|) with $ replaced by $ s . 

(ii) Evolution of $ a is gwen by the equation $ a = L4,$ a ]. 

(iii) If n = l, then $ a = 0. 

Proof. The matrix M l satisfies the same equation Q as does M 

M* = JH"(x)M t + M t H"(x)J t . 

Therefore © is satisfied separately by M s and M a . Correspondingly, the second equation 
in (fTUj) is satisfied separately by $ s and $ a . 

The following calculation shows that the RHS of the first equation in (jl(lj) depends only 
on $ s : 

tr = tr (iJ"M) = tr (^"M s ) = tr 

Step (*) holds because H" is symmetric and tv(H"M) = tr(M*if") = tr(i?"M*). Finally, 
in case n = 1 the matrix $ a is proportional to I2 , so $ a = [A, I2] = 0. □ 

We are interested in an influence of the described extension on the original dynamics 
due to the entanglement of the x and $ variables. From this point of view, it suffices to 
consider system (jlO)) assuming that matrix $ has the form $ = JM with a symmetric M. 
Consequently, it obeys $'J+ J$ = 0. Recall that In x 2n real matrices with this property 
are called symplectic. They form the classical Lie algebra sp(2n, M). 

Definition 3. System (jlUj) with $ £ sp(2n, M.) is called the two-system of Hamilton's 
system (jTJ). 

Remark. The term two-system and a family of n-systems were introduced in |BK 19 94]. 
According to |BK 1994j . n-system refines the classical limit approximation for a localized 
solution of the Schrodinger equation by taking into account all moments of the wave packet 
up to order n. Higher n-systems seem not to be practical because they add too many 
degrees of freedom to the original system. A Hamiltonian formulation has not been found 
for n-systems with n > 2. The situation is quite different for the 2-system, as we discuss in 
the next section. 
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4. Hamiltonian structure of the two-system 

a. By definition, the vector form (J2J) is a Hamiltonian system. However, it is not immediately 
obvious which Hamiltonian structure (if any) is related to to the two-system. Such a structure 
was found in [S 1994j. The result will be recalled in d; it uses the notion of a linear Poisson 
bracket on sp(2n, R) (more precisely, on its dual space, but we don't make a difference using 
the standard inner product on the space of matrices: (A, B) = tr(AB*) ). 

Preliminary, in b, c we describe another Hamiltonian approach to the two-system, con- 
sidering a family of multivector forms, which generalize the vector form. The phase space 
of the two-system is decomposed into invariant subsets corresponding to fixed ranks and 
signatures of the symmetric matrix M = — J<3>. 

b. Let m + , m_ and m be the numbers of positive, negative, and zero eigenvalues of a 
symmetric 2n x 2n matrix M. Then 

m + + m_ + m = 2n, m + + m_ = rk M. 

The pair (m + ,m_) is called signature of M. 

Proposition 4. Let a symmetric matrix-function M(t) satisfy the equation 

JM = [JM, A{t)\ 

with a function A : R — > sp(2n, R). Then signature of M(t) does not change. 

Proof. There exists a non-degenerate matrix function B(t) such that 

JM(t) = B{t) JM(0) Bit)- 1 

Multiplication by a non-degenerate matrix preserves rank, so rkM(t) = const. The ma- 
trix function JM(t) is continuous and it has constant rank and real eigenvalues. Nonzero 
eigenvalues can not turn to zero, hence their signs don't change. □ 

Definition 4. Define a subset X(m + ,m_) of the phase space M 2n x sp(2n, R) as the 
set of all pairs (x, JM) with matrix M of signature (m + ,m_). It is invariant (Proposition 
4). Restriction of the two-system on X(m + ,m_) is called the two-system in signature 
(m + , m_) . 

c. Construction (cf. Sect.EJ). Let m + , m_ be nonnegative integers and r = m + + m_ < 2n. 
Take the sum of r + 1 copies of M 2n to be the phase space of a new dynamical system. Let x, 
{y;} — yij • • • ^ y m+) { z j} = zi, ■ ■ ■ , z m_ be the respective groups of dynamical variables. 
Equip the phase space with the standard symplectic structure J © J © • • • © J and define a 
new Hamilton's function 

W(x, {y,}, {z,}) = H(x) + £F(x, y 4 ) - ]TF(x, Zj ), (12) 

i=i j=i 

where H(x) is Hamilton's function of system (0) and F is defined in 
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Definition 5. The system described in the Construction is called the multi-vector form 
of signature (m + ,m_) (for a Hamiltonian combination of the systems (JIJ and ©))• For 
brevity, we'll write multivector (m+,m-)-form. 

A symmetric matrix M of signature (m + ,m_) can be written in the form 

m+ rn- 

M = $>, yJ - 5>z* (13) 

3=1 i=l 

The matrix built from the vector-columns y » / 1 1 y » 1 1 , z j/|| z j|| is determined up to a pseudo- 
orthogonal rotation, namely to an element of the group 0(m + ,m_) x I mo . Consequently, 
the two-system in signature (m + , m_) can be obtained from the multivector (m + , m_)-form 
by projection. A precise formulation follows. 

Proposition 5. Let (x , JM ) be the initial point of a trajectory (x(i), JM(t) ) of the 
two-system in the signature (m + ,m_). Let {yi(0)}™ + ; {zj(0)}™~ be the components of 
some decomposition of M of the form (jrHj) . If (x(t) ; {y«(t)} ; {zj(t)}) is the solution of 
the multivector {m + ,m J) -form with initial data (xo, {yi(0)} ; {zj(0)}) ; then for all t 

x(t) = x(t) 

and 

M(t) = Eyi(*)y5(*) - E 

i=l 3=1 

Proof. Assuming M(t) in the said form, it obeys if Vj and z» obey (j2J). Also the evolution 
of x(t) determined by the Hamiltonian (|12|) coincides with that in the first equation (jlOj) . 

□ 

Note that a multivector form cannot be subdivided into vector forms (J7J) since x(t) is 
subject to a collective effect by all y's and z's. 

d. The two-system also possesses a non-trivial Hamiltonian structure with a degenerate Pois- 
son bracket. (Note without going into details that this fact reflects some general construction 
of symplectic geometry |KM 1 996, § 1.5].) 

Identify the phase space of the two-system with the space R 2n x sp(2n, M). For functions 
on this space, introduce a Poisson bracket by the formula 

{f/,V}(x,$) = U'JV' + 2tr(V [V$C/,V$y]) . (14) 

Here U' denotes the gradient of U(x, $) with respect to x-variables, and V$Z7 denotes the 
matrix from sp(2n, R) with entries 

Oil 

In particular, 

V$(det$) = ($- x )*det$; V»(tr($A) ) = A\ VA (15) 

The bilinear operator (|14|) is indeed a Poisson bracket, i.e. it is anti-symmetric and satisfies 
the Leibnitz and Jacobi identities. 
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Proposition 6. The two-system ([10)1 can be written in the Hamiltonian form 

Xi = {x h W(x, $)}, $y = , W(x, $)} 
wzi/i Poisson bracket fj!4|) and Hamilton's function 

W(x; $) = ff(x) - itr(i4(x)$), A(x) = H"(x)J 



(16) 



(17) 



Proof. Equivalence of the first equations in ()16j) and (fTUj) immediately follows from ()14j) . 
Then, V$$ij = i£y , the matrix with the only nonzero (unity) entry [ij]. The second 
equation in (JTSJ) gives = -(1/2) A 1 . Now (JUJ) yields 



□ 



W(x, $)} = tr (V + A%)J = tr(% [A, $]*) = [A, 

so the right-hand sides of the second equations in (fTt^) and (JTUJ) agree. 



Remark. The Poisson bracket (J 14)) can be described via the matrix function Q(x, $) on 
the (2n 2 + 3n)-dimensional linear space R n x sp(2ra, R). The matrix at a generic point 
has corank n. Moreover, we have the following 

Proposition 7. Besides Hamilton's function 7i, the two-system has n additional, func- 
tionally independent integrals, which are Casimir's functions of the bracket (|14jl. As such 
integrals, one can take the coefficiens of the characteristic polynomial P(\) = det($ — XI) 
at A , A 2 , . . . , A 2(n_1) . 

Proof. Note that det $ is in involution with any function Z7(x, $) with respect to the bracket 

{det$, ^7(x,$)} = 0, 

which can be shown using the first of equations (|15|). 

Since $ is symplectic, its characteristic polynomial is an even function of degree 2n. So 
it has n non-trivial coefficients corresponding to the powers 0, 2, ... , 2(n — 1) of A. They 
can independently take any values. □ 

It is known from the general theory of Hamiltonian systems with degenerate Poisson 
brackets (see e.g. |G 2001j . |DS 1989] ) that the phase space of the 2-system is a disjoint 
union of invariant subspaces — symplectic leafs. Restriction of the bracket on any symplectic 
leaf is non-degenerate, hence, at least locally, there exist canonical coordinates on the leaf, 
in which the system has the standard Hamiltonian form. A point of the phase space in a 
generic posisiton belongs to a symplectic leaf of the lowest possible codimension n and the 
corresponding nondegenerate Hamiltonian system has n 2 + n degrees of freedom. There 
exist symplectic leafs of higher codimensions (but always of an even dimension). 

Example: Case n = 1. The two-system has order 5: two original variables q, p, and three 
"moments of the second order" a, (3, 7, cf. the system in the Introduction. The explicit 
expression for the matrix Q is given by 



VI = J © u, J 






1 " 











-1 



2a Af3 
-2a 27 
-4/3 -2 7 
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The (degenerate) Poisson bracket on the phase space is {U, V} = J2i.j=i The 
bracket has one Casimir function 

S = a/3-7 2 . (19) 
5. Integrable cases and special solutions 

Integrability of the two-system (fTUj) has a rare occurence. We discuss two cases where it 
takes place. 

Proposition 8. Two-system of a Hamiltonian system with Hamilton's function being a 
polynomial of degree 2 or less splits into two independent subsystems and has the general 
solution in a closed form. 

Proof. For any i,j, k, H"- k = 0, and the matrix H" of the second derivatives is constant. 
Thus the system ((TP)) splits into the original system and the system for <3>(i) with constant 
matrix A. Both subsystems are linear, have constant coefficients, and therefore integrable 
in elementary functions. □ 

Proposition 9. Two-system of an integrable Hamiltonian system written in action-angle 
variables, with Hamilton's function depending only on the action variables, is integrable. 

Proof. We'll consider only case n = 1, where the nature of integrability is same as in the 
general case. Let (/, 9) be the action-angle variables, and H = H(I) Hamilton's function. 

Denote M = ( ^ ^ V The two-system reads 

H'(I) + H"'(I)a, a = 0, p = H"{I)p, 7 = 2H"(I)f3 (20) 

The general solution of (J20|) is easily found to be 

I = I , a = a , = O o + tui, P = Poe UJ2t , 7 = Wo{e W2t - 1) + 7o- 

Here Jo, #0, «o, Po, 7o are arbitrary constants and lo\ = H'(I Q ) + H"'(I )ao, uj 2 = H"(I ). 

□ 

Proposition 9 shows that the two-system in action-angle variables is always integrable, 
while in general it is non-integrable except in the situation of Proposition 8. Thus the two- 
system is essentialy not invariant or covariant under canonical transformations of system (JTj). 
Neveretheless, there are two cases where the x-component of a solution of the two-system 
obeys equation (JTJ). 

Proposition 10. (i) Consider special initial data for the two-system 

x(0) = x , $(0) = 0. 



a 


P 


P 


7 


0. 
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The solution of (|10|1 with these initial data is (x(i|xo), $ = 0), where x(t|xo) is the solution 
of (JTJ) with initial data x . 

(ii) Let Xq be such a stationary point of system (QP where all third derivatives of the Hamil- 
tonian vanish: H"L{xq) = 0. Then the two-system (fTU|) has a solution, for which x(t) = x Q 
and obeys the linear system with constant coefficients 

$=[A(x ),$\. (21) 
Proof easily follows from the structure of system fllOl) . □ 
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